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Abstract. We study the kinetic and chemical equihbration in 'infinite' parton-hadron matter within the Parton- 
Hadron-String Dynamics transport approach, which is based on a dynamical quasiparticle model for partons 
matched to reproduce lattice-QCD results - including the partonic equation of state - in thermodynamic equi- 
librium. The 'infinite' matter is simulated within a cubic box with periodic boundary conditions initialized at 
different baryon density (or chemical potential) and energy density. The transition from initially pure partonic 
matter to hadronic degrees of freedom (or vice versa) occurs dynamically by interactions. Different thermody- 
namical distributions of the strongly-interacting quark-gluon plasma (sQGP) are addressed and discussed. 



1 Introduction 



2 The model 
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Nucleus-nucleus collisions at ultra-relativistic energies are 
studied experimentally and theoretically to obtain informa- 
tion about the properties of hadrons at high density and/or 
temperature as well as about the phase transition to a new 
state of matter, the quark-gluon plasma (QGP). Whereas 
the early 'big-bang' of the universe most likely evolved 
through steps of kinetic and chemical equilibrium, the lab- 
oratory 'tiny bangs' proceed through phase-space configu- 
rations that initially are far from an equilibrium phase and 
then evolve by fast expansion. On the other hand, many 
observables from strongly-interacting systems are domi- 
nated by many-body phase space such that spectra and 
abundances look 'thermal'. It is thus tempting to charac- 
terize the experimental observables by global thermody- 
namical quantities like 'temperature', chemical potentials 
or entropy |1 23456 78]. We note, that the use of macro- 
scopic models like hydrodynamics |9 10 11 12 1 employs 
as basic assumption the concept of local thermal and chem- 
ical equilibrium. The crucial question, however, how and 
on what timescales a global thermodynamic equilibrium 
can be achieved, is presently a matter of debate. Thus nonequi- 
librium approaches have been used in the past to address 
the problem of timescales associated to global or local equi- 
libration 113 14 15 16 17 18 19 20 21 1. In view of the in- 
creasing 'popularity' of thermodynamic analyses a thor- 
ough microscopic study of the questions of thermalization 
and equilibration of confined and deconfined matter within 
a transport approach appears necessary. 
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In this contribution, we study the kinetic and chemical equi- 
libration in 'infinite' parton-hadron matter within the novel 
Parton-Hadron- String Dynamics (PHSD) transport approach 
r22'231, which is based on generalized transport equations 
on the basis of the off-shell Kadanoff-Baym equations 124 25)1 
for Green's functions in phase-space representation (in the 
first order gradient expansion, beyond the quasiparticle ap- 
proximation). In the KB theory, the field quanta are de- 
scribed in terms of propagators with complex self-energies. 
Whereas the real part of the self-energies can be related 
to mean-field potentials, the imaginary parts provide in- 
formation about the lifetime and/or reaction rates of time- 
like "particles". The basis of the partonic phase descrip- 
tion is the dynamical quasiparticle model (DQPM) 12612711 
matched to reproduce lattice QCD results - including the 
partonic equation of state - in thermodynamic equilibrium 
|24|. In fact, the DQPM allows a simple and transparent 
interpretation for thermodynamic quantities as well as cor- 
relators - measured on the lattice - by means of effec- 
tive strongly interacting partonic quasiparticles with broad 
spectral functions. The transition from partonic to hadronic 
degrees of freedom is described by covariant transition rates 
for fusion of quark-antiquark pairs or three quarks (anti- 
quarks), obeying flavor current conservation, color neutral- 
ity as well as energy-momentum conservation. 

The 'infinite' matter is simulated within a cubic box 
with periodic boundary conditions initialized at various val- 
ues for baryon density (or chemical potential) and energy 
density. The size of the box is fixed to 9^ fm^. The initial- 
ization is done by populating the box with light (u, d, s) 
quarks, antiquarks and gluons with random space posi- 
tions and the momenta distributed according to the Fermi- 
Dirac distribution. The total numbers of the quarks and an- 



EPJ Web of Conferences 



hadrons 
partons 



> 
O 

— 

z 

■a 





Fig. 1. Snapshot of the spatial distribution of hadrons (blue) and 
partons (red) at an evolution time of 40 fm/c after the systems was 
initialized by solely partons at an energy density below critical. 
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Fig. 3. The reaction rates for elastic parton scattering (blue), 
gluon splitting (green) and flavor neutral qq fusion (red) as a func- 
tion of time. 




tiquarks are chosen so that the system with various desired 
values of the energy density s (the total energy of the par- 
ticles divided by the size of the box) and baryon potential 
fi can be studied. 

In the course of the subsequent transport evolution of 
the system by PHSD, the numbers of gluons, quarks and 
antiquarks are adjusted dynamically through the inelastic 
collisions to equilibrium values, while the elastic collisions 
lead to eventual thermalization of all the particle species 
(e.g. u, d, s quarks and gluons, if the energy density in the 
system is above critical). Please note that if the energy den- 
sity in a local cell drops below critical either due to the lo- 
cal fluctuations or because the system was initialized with 
a low enough number of partons, a transition from initial 
pure partonic matter to hadronic degrees of freedom occurs 
dynamically by interactions. 



3 Results 

We present in the Fig. [T] a snapshot of the spatial distri- 
bution of hadrons (blue) and partons (red) at an evolution 
time of 40 fm/c after the systems was initialized by solely 
partons at an energy density of 0.4 GeV/fm^. At this en- 
ergy density - slightly below the critical energy density for 
the deconfinement phase transition - most of the partons 
have already formed hadrons at this point of the the evolu- 
tion of the system, but a small fraction of partons have yet 
not hadronized. The remaining partons are approximately 
in the thermal equilibrium with the hadrons. 

In Fig. |2] we show a snapshot of the systems that has 
been initialized at an energy density of 2. 1 8 GeV/fm^, which 
is clearly above the critical energy density. One can see 
in Fig. |2] the spatial distribution of light quarks and an- 
tiquarks (red), strange quarks and antiquarks (green) and 
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Fig. 5. Abundances of the u,d s quarks + an tiquarks and gluons as 
a function of time for a system initialized at an energy density of 
9.43 GeV/fm^. 
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Fig. 6. The average abundances of hadrons (blue) and partons 
(red) in equilibrium as functions of the energy density e. 



gluons (blue) at a time of 40 fm/c. At this energy density 
no hadrons are seen. The different parton species are in 
thermal and chemical equilibrium. 

After a few fm/c the system - initialized at an energy 
density of 2.87 GeV/fm^ - has achieved chemical and ther- 
mal equilibrium, since the reactions rates are practically 
constant and obey detailed balance for gluon splitting and 
qq fusion. This is shown in Fig.[3l where the reaction rates 
for elastic parton scattering (blue), gluon splitting (green) 
and flavor neutral qq fusion (red) are presented as a func- 
tion of time. 

Another indication that the system has achieved the 
thermal equilibrium is seen in the distribution of the ki- 
netic energy of the particles. We show in Fig. |4] the en- 
ergy spectra for the off-shell u (red) and s quarks (green) 
and gluons (blue) in equilibrium for a system initialized 
at an energy density of 5.37 GeV/fm^. The spectra may 
well be described by a Boltzmann distribution with tem- 



perature T=240 MeV in the high energy regime. The devi- 
ations from the Boltzmann distribution at low energy E are 
due to the broad spectral functions of the partons. 

On the other hand, a sign of the chemical equilibra- 
tion is the stabilization of the abundances of the different 
species. In Fig. |5] we show the particle abundances as a 
function of time for a system initialized at 9.43 GeV/fnr'. 
One can see that the chemical equilibration is reached after 
about 15 fm/c. 

It is interesting to observe in Fig. |6] the average abun- 
dances of hadrons (blue) and partons (red) in equilibrium 
as functions of the energy density s. At energy densities 
below critical 0.5 GeV/fm-') the system evolves into a 
state, which is dominated by hadrons and has a very small 
fraction of partons due to rare fluctuations of local energy 
density to high values. At higher energy densities, the sys- 
tem is in a QGP final state, with a small hadron admixture. 
At high enough energy (above approx. 2 Ge V/fm^*) we find 
that hadron fraction is negligible. In the regime of energy 
densities from 0.48 to 1.3 GeV/fm-' the calculations have 
provided so far no stable equilibrium over time due to large 
fluctuations between hadronic and partonic configurations. 
Further studies are on the way. 



4 Conclusions 

We have studied the kinetic and chemical equilibration in 
'infinite' parton-hadron matter within the Parton-Hadron- 
String Dynamics transport approach (PHSD), which is based 
on a dynamical quasiparticle model for partons (DQPM) 
matched to reproduce lattice-QCD results - including the 
partonic equation of state - in thermodynamic equilibrium. 

The 'infinite' matter has been simulated within a cubic 
box with periodic boundary conditions initialized at dif- 
ferent baryon density (or chemical potential) and energy 
density. 

Depending on the energy density, the system evolved 
into an ensemble of either hadrons or partons in chemical 
and thermal equilibrium. Abundances of the final particles 
depend on the energy density. The temperature of the de- 
grees of freedom in the final state was measured by fitting 
the slopes of the Boltzmann-like tails of the thermal distri- 
butions of their kinetic energy. 
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